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We propose a novel spin liquid state for a spin 5 = 1 antiferromagnet in two dimensions. The 
ground state violates P and T, is a spin-singlet, and is fully invariant under the lattice symmetries. 
The spinon and holon excitations are deconfined and obey non-abelian statistics. We present pre- 
liminary numerical evidence that the universality class of this topological liquid can be stabilized 
by a local Hamiltonian involving three-spin interactions. We conjecture that spinons in spin liquids 
with spin larger than 1/2 obey non-abelian statistics in general. 
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Introduction — Fractional quantization in two-dimen- 
sional quantum liquids is witnessing a renaissance of in- 
terest in present times. The field started about a quar- 
ter of a century ago with the discovery of the fractional 
quantum Hall effect, which was explained by Laughlin [lj 
in terms of an incompressible quantum liquid support- 
ing fractionally charged (vortex or) quasiparticle excita- 
tions. When formulating a hierarchy of quantized Hall 
states 0, IH to explain the observation of quantized Hall 
states at other filling fractions fractions, Halpcrin J3[ 
noted that these excitations obey fractional statistics [4| , 
and are hence conceptually similar to the charge-flux 
tube composites introduced by Wilczek two years ear- 
lier [|. 

The interest was renewed a few years later, when An- 
derson proposed that hole-doped Mott insulators, and 
in particular those described by the t-J model universally 
believed to describe the CuO planes in high T c supercon- 
ductors, can be described in terms of a spin liquid (i.e., 
a state with strong, local antiferromagnetic correlations 
but without long range order), which would likewise sup- 
port fractionally quantized excitations. In this proposal, 
the excitations are spinons and holons, which carry spin 
1/2 and no charge or no spin and charge +e, respectively. 
The fractional quantum number of the spinon is the spin, 
which is half integer while the Hilbcrt space (for the un- 
doped system) is built up of spin flips, which carry spin 
one. One of the earliest proposals for a spin liquid sup- 
porting deconfined spinon and holon excitations is the 
(abclian) chiral spin liquid (CSL). Following up on an 
idea by D.H. Lee, Kalmcycr and Laughlin |7[ proposed 
that a quantized Hall wave function for bosons could be 
used to describe the amplitudes for spin-flips on a lattice. 
The CSL state did not turn out to be relevant to CuO 
superconductivity, but remains one of very few examples 
of two-dimensional spin liquids with fractional quantiza- 
tion. The other established examples are the resonat- 
ing valence bond (RVB) phases of the Rokhsar-Kivelson 
model Q on the triangular lattice identified by Moessner 
and Sondhi and of the Kitaev model [13] ■ 

The present renaissance of interest in fractional quan- 
tization is due to possible applications of states support- 



ing excitations with non-abelian statistics to the rapidly 
evolving field of quantum computation and cryptogra- 
phy. The paradigm for this universality class is the Pfaf- 
fian state introduced by Moore and Read ll| in 1991. 



The state was proposed to be realized at the experimen- 
tally observed fraction v = 5/2 [l2[ [i.e., at v = 1/2 in 
the second Landau level) by Wen, Wilczek, and one of 
us fl3| . a proposal which recently received experimental 
support through the direct measurement of the quasipar- 
ticle charge [TJ, [l5|. Pfaffian type states are further con- 
jectured to be realized for one-dimensional bosons with 
three-body hard core interactions in general [lij ]. The 
Moore-Read state possesses p + ip wave pairing correla- 
tions. The flux quantum of the vortices is one half of the 
Dirac quantum, which implies a quasiparticle charge of 
e/4. Like the vortices in a p wave superfluid, these quasi- 
particles possess Majorana-fermion states [17[ at zero en- 
ergy (i.e., one fermion state per pair of vortices, which 
can be occupied or unoccupied). A Pfaffian state with 2L 
spatially separated quasiparticle excitations is hence 2 L 
fold degenerate, in accordance with the dimension of the 
internal space spanned by the zero energy states. While 
adiabatic interchanges of quasiparticles yield only over- 
all phases in abelian quantized Hall states, braiding of 
half vortices of the Pfaffian state will in general yield 
non-trivial changes in the occupations of the zero en- 
ergy states IH lH, which render the interchanges non- 
commutative or non-abelian. In particular, the inter- 
nal state vector is insensitive to local perturbations — it 
can only be manipulated through braiding of the vor- 
tices. These properties together render non-abelions pre- 
eminently suited for applications as protected qubits in 
quantum computation [2fJ. Non-abelian anyons further 
appear in certain other quantum Hall states including the 
Read-Rezayi states [2l[ , in the Kitaev model [l(J , and in 
the Yao-Kivelson model [22| ■ 

In this Letter, we propose a novel chiral spin liquid 
state for an S = 1 antiferromagnet. The spinon and 
holon excitations of this state are deconfined and obey 
non-abelian statistics, with the braiding governed by Ma- 
jorana fermion states. The state violates time reversal 
(T) and parity (P), is a spin singlet, can be formulated on 



any lattice type, and fully respects all the lattice symme- 
tries. The state possesses a 3-fold topological degeneracy 
on the torus geometry. We provide preliminary numeri- 
cal evidence that the state can be stabilized on the trian- 
gular lattice by a local Hamiltonian involving three-spin 
interactions. Finally, we conjecture that spinons in spin 
liquids with spin larger than 1/2 might obey non-abelian 
statistics in general. 

Non-abelian chiral spin liquid state — The state we pro- 
pose is most easily written down for a circular droplet 
with open boundary conditions occupying N sites of a 
triangular or square lattice 5 = 1 antiferromagnct. The 
wave function for re-normalized spin flips, 
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is given by a bosonic Pfaffian state in the complex coor- 
dinates z = x + iy supplemented by a gauge phase G(ij a ). 
The Pfaffian is given by the fully antisymmetrized sum 
over all possible pairings of the N coordinates, 
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The "particles" Zi represent re-normalized spin flips act- 
ing on a vacuum with all spins in the S z — — 1 state, 
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where the sum extends over all possibilities of distribut- 
ing the N "particles" over the N lattice sites allowing for 
double occupation, and 
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The lattice may be anisotropic; we have chosen the lattice 
constants such that the area of the unit cell spanned by 
the primitive lattice vectors is set to unity. For a triangu- 
lar or square lattice with lattice positions given by rj n ^ m = 
na + mb, where a and b are the primitive lattice vectors 
in the complex plane, G{r) n<m ) = 

Singlet property — While the topological properties, 
and in particular the non-abelian statistics of the frac- 
tionalized excitations of ([1]), are suggestive to those fa- 
miliar with Pfaffian states, the invariance under spin ro- 
tation and lattice symmetries is less so. We content our- 
selves here with a direct proof of the singlet property, 
which at the same time serves to motivate the necessity 
for the re-normalization of the spin- flip operators ([4]). 

Since S* ot \ipo) = by construction, it is sufficient to 
show Sfo t \ipo) = 0. Note first that when we substitute ((T|) 
with ^ into we may omit the antisymmetrization 
A in as it is taken care by the commutativity of the 
bosonic operators S a . (Throughout this Letter, we do 



not keep track of overall normalization factors.) Let ipo 
be ipo without the operator A in ((2|). Since ipo[ z i] is still 
symmetric under interchange of pairs, we may assume 
that a spin flip operator S~ acting on |^o) will act on 
the pair (zi, z<z): 
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where we have used 
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This implies 5 t ~ t |V>o) = Z) a =i $a IV'o) = if and only if 
Ea=i i>(Voi, Z2 1 za, ■ ■ •) = V z 2 ,zs, ■ ■ -z N . The Perelo- 
mov identity [24j states that this holds for lattice sums 
of e~5' ,? °' G(r] a ) times any analytic function of r] a . 

Generation from filled landau levels — Rather than pro- 
ceeding in verifying invariance properties of the non- 
abelian CSL state (H|) directly, wc motivate them indi- 
rectly through demonstrating that the state can alterna- 
tively be generated though successive projection via the 
abclian CSL from the wave functions of a filled lowest 
Landau level (LLL). If wc choose an auxiliary magnetic 
field with a strength of one half of a Dirac flux quanta 
per lattice site, the wave function for a circular droplet 
of M = £■ fermions filling the LLL is given by 
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The (abclian) CSL state for spin S = | [tJ, which was 
recently shown to be the unique and exact ground state 
of a local Hamiltonian [2511 , 
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where the "particles" Z{ describe spin flips 5+ acting 
on a "vacuum" state with all the spins |, and G(r/ a ) 
is as above, can be generated by Gutzwiller projection 
of the LLL |5]) filled once with f and once with J, spin 
fermions 1261. 12711 : 
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where the sum extends over all partitions of the lattice 
sites into z's and w's and the cAs are fermion creation 
operators. We can rewrite the CSL state vector in terms 
of Schwinger bosons A and , 
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provided we define \to SL [cj, cj] sucn that tne operators 
are ordered according to a fixed labeling of the lattice 
sites. The non-abelian CSL state (Q]) can thus alterna- 
tively be written as a symmetrization over two abelian 
CSL states 
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where 5 indicates symmetrization. Since the LLL states 
([5]) are (on compact surfaces) translationally and rota- 
tionally invariant modulo gauge transformations in the 
auxiliary magnetic field, and ([7]) is manifestly gauge co- 
variant, both the CSL states © and (Q]) are invariant 
under lattice transformations. Note that this projective 
construction also implies the singlet property of the CSL 
states. It can be used to formulate the CSL states on any 
lattice, and to generalize them to arbitrary spin: 
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Written in terms of (then differently) re-normalized spin 
flip "particles" , the wave function generalizes from a 
bosonic Pfaffian state for S = 1 to bosonic Read-Rczayi 
states for S > 1. 

Non-abelian spinon and holon excitations — The spinon 
excitations of (Q]) are analogous to the half vortex quasi- 
particles of the Moore- Read quantum Hall state [ll| . For 
example, to create 4 J. spin spinons at locations r}\, 772, 773, 
and r)i, we simply insert half quantum vortices inside the 
Pfaffian @, which then becomes 
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The braiding properties of the spinons are insensitive to 
the spinon spin, and are exactly those of the Moore- Read 
quasiparticles [o, 0, • The proof of the singlet prop- 
erty given above can be extended to show that a pair of 
I spin spinons transforms as an S = 1 triplet excitation, 
which implies that each spinon carries spin S = With 
the implicit assumption that the S = 1 spins on each 
lattice site consist of two electrons in triplet configura- 
tions, we can create holon excitations by annihilating I 
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FIG. 1: (Color online) The eleven interaction terms included 
in our trial Hamiltonian with the numerically optimized co- 
efficients (see text). 

spin electrons on sites with J, spin spinons. The braiding 
properties of the holons are equivalent to those of the 
spinons. 

Model Hamiltonian — The first question with regard to 
possible applications of our state to quantum compu- 
tation is whether a state belonging to the universality 
class described by ((l) can be stabilized through a local 
Hamiltonian. While we are short of a definite answer, 
we have done our best to address the question numeri- 
cally. To begin with, we have written out the state (JXJ) for 
an isotropic, triangular lattice with 16 sites and periodic 
boundary conditions, which imply a three- fold topologi- 
cal degeneracy [l3|. We then numerically optimized the 
coefficients of a set of local spin interaction terms (see 
Fig. [I} such that the ground state of our trial Hamil- 
tonian is energetically closest to a suitable linear com- 
bination of the three (in the thermodynamic limit de- 
generate) Pfaffian states, which we then compare to the 
exact eigenstates. As shown in Fig. [21 the three lowest 
energy eigenstates of our trial Hamiltonian have a signif- 
icant overlap with the Pfaffian states (i.e., 0.959, 0.964, 
and 0.934 in a fully symmetry reduced 5f ot = Hilbcrt 
space with dimension 163101), which suggests that the 
exact states belong to the same universality class. Note 
that the coefficients in Fig. Q] fall off rapidly with the 
distance. Small variation of the parameters induce no 
sensitive change in the overlaps, which indicates that the 
non-abelian CSL state is stablilized throughout a finite 
region in parameter space. Our evidence is unfortunately 
not conclusive as the three CSL states are not separated 
by a large gap from the remainder of the spectrum, which 
indicates that the system we can access numerically is too 
small to settle the question unambiguously. 

Experimental realization — Recent work on polar 
molecules in optical lattices p8|, b ut in particular on 
engineering 3-body interactions 29j, suggests that a re- 



3 




overlaps 



|*0> 



H 1 1 1 h- 

r A X K M 



-+- 
s 



-0.90 
-0.92 
-0.94 
-0.96 
-0.98 
-1.00 
-1.02 
-1.04 



FIG. 2: (Color online) Spectral plot of our trial Hamiltonian 
in comparison with the energy expectations values for the 
three (in the infinite system topologically degenerate) Pfaffian 
ground states states at the F point. The inset shows the 
overlap of the Pfaffian states with the three lowest states of 
our Hamiltonian. 



the Hilbcrt space spanned by the ground state plus all 
states with different numbers of spinons for the spin liq- 
uid we propose is 3^ , as required for a 5 = 1 system with 
N sites. We conjecture that the non-ablelian statistics 
for the spinons is not only a sufficient, but even neces- 
sary condition for the state counting to work out consis- 
tently. (Haldane [32| has shown that the state counting 



for 5 



spin liquids works out consistently if one as- 



sumes abelian half-fcrmi statistics for the spinons.) 

Conclusion — In this work, we have constructed an 
5 = 1 CSL and argued that its spinon and holon exci- 
tations obey non-abelian statistics. We have used exact 
diagonalization studies to obtain preliminary indication 
that the state can be stabilized on a 5 = 1 triangular 
lattice. 

MG wishes to thank B. Paredes for a highly stimulating 
discussion. RT was supported by a PhD scholarship from 
the Studienstiftung des deutschen Volkes. 



alization of the non-abelian CSL proposed here might be 
possible at some stage in the future. 

Non-abelian spinons in general — Efforts to understand 
high T c superconductivity in terms of an RVB spin liq- 
uid have revealed a general connection between d-wave 
superconductors and 5=4 spin liquids on the square 
lattice HSU. In particular, a wide class of (undoped) 



5 = h spin liquids can be obtained by Gutzwiller projec- 
tion from the wave function of a <i-wave superconductor 
with suitably chosen parameters. This suggests a general 
connection between the (abelian) vortices of the super- 
conductor and the (abelian) spinons in the spin liquid. 
If one Gutzwiller projects a d + id wave superconductor 
with suitably chosen parameter on a square lattice, one 
obtains exactly the CSL state ((6J. 

The p + ip pairing correlations in the non-abelian CSL 
state (HJ) introduced above suggest a similar correspon- 
dence between the non-abelian vortices of the supercon- 
ductor and the non-abelian spinon excitations (|10p . As 
in the abelian case 5 = \ , the P and T violation of the 
state appears to be necessary for the spinon to be de- 
confined, but does not seem essential to the topological 
properties. We are hence led to conjecture that there is a 
general connection between p-wave superfluids and 5 = 1 
spin liquids, in that the non-abelian braiding properties 
of the vortices of the superfluid are also general prop- 
erties of the spinons in 5 = 1 antiferromagnets. True, 
the spinons will only be free under special circumstances, 
and the propensity to be confined will only increase with 
the spin 5. Even in an ordered antiferromagnet, how- 
ever, spinons (and holons) are the fields appropriate for 
describing the physics at sufficiently high energy scales, 
i.e., energies above the ordering temperature. 

We will show elsewhere that the total dimension of 
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